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Diffractive Higgs boson production at the Tevatron and LHC 

R. Enberg, 1 G. Ingelman, 1,2 A. Kissavos, 1 and N. Timneanu 1 

1 High Energy Physics, Uppsala University, Box 535, S-75121 Uppsala, Sweden 
2 Deutsches Elektronen- Synchrotron DESY, Notkestrasse 85, D-22603 Hamburg, Germany 

Improved possibilities to find the Higgs boson in diffractive events, having iess hadronic activity, 
depend on whether the cross section is large enough. Based on the soft color interaction models 
that successfully describe diffractive hard scattering at HERA and the Tevatron, we find that only 
a few diffractive Higgs events may be produced at the Tevatron, but we predict a substantial rate 
at the LHC. 
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The Higgs boson is predicted as the physical manifesta- 
tion of the mechanism giving masses to the fundamental 
particles in the Standard Model. The discovery of this 
missing link is of top priority in particle physics. Based 
on the discovery ||, |] of diffractive hard scattering pro- 
cesses H it has been considered whether the Higgs can be 
more easily observed in diffractive events at high energy 
hadron colliders. The lower hadronic activity in such 
events with large rapidity gaps should improve the possi- 
bilities to reconstruct the Higgs from its decay products. 

In single diffraction a beam hadron emerges quasi- 
elastically scattered with a large fraction of the orig- 
inal momentum and separated by a gap in polar an- 
gle (pseudorapidity 77 = —In tan 0/2) from a produced 
hadronic system X. Particularly clean events are pro- 
duced in so-called double pomeron exchange (DPE), 
where both beam hadrons emerge intact separated by 
rapidity gaps from a centrally produced X-system. An 
extreme possibility is exclusive Higgs production, pp — > 
ppH, where the central system is just a Higgs boson that 
may be reconstructed using a missing mass method [Q. 
The crucial question for the usefulness of these diffractive 
Higgs production processes is whether their cross sections 
are large enough. 

Predicted cross sections vary by orders of magnitude 
between calculations based on different models || . Some 
predictions for the Fermilab Tevatron are large enough 
to be of experimental interest, whereas others are not. 
The more limited energy, compared to LHC, implies a 
stronger kinematical suppression to produce the heavy 
Higgs boson in such an A-system, having only a frac- 
tion of the overall invariant mass of the collision. In this 
Letter we improve on this theoretical uncertainty by pre- 
senting results on diffractive Higgs production based on 
the recently developed soft color interaction models. In 
contrast to other models used for estimating the diffrac- 
tive Higgs cross section, our models have proven very 
successful in reproducing experimental data on diffrac- 
tive hard scattering processes both from the DESY ep 
collider HERA and from pp collisions at the Tevatron 
||. This puts us in a good position to give predictions 
on diffractive Higgs production. 

The soft color interaction (SCI) model and the 



generalized area law (GAL) model || were developed 
in an attempt to better understand non-perturbative 
QCD dynamics and provide a unified description of all 
final states. The basic assumption is that soft color ex- 
changes give variations in the topology of the confining 
color string-fields which then hadronizc into different fi- 
nal states, e.g. with and without rapidity gaps or leading 
protons. Also other kinds of experimental results are 
described in a very economical way with only one new 
parameter. Particularly noteworthy is the turning of a 
cc pair in a color octet state into a singlet state producing 
charmonium H in good agreement with observed rates. 




FIG. I: Higgs production in pp collisions with string topolo- 
gies (double-dashed lines) before and after soft color interac- 
tions in the SCI or GAL model, resulting in events with one 
or two rapidity gaps (leading particles). 

The SCI model Q is implemented in the Lund Monte 
Carlo programs lepto for deep inelastic scattering 
and pythia [[ll] for hadron-hadron collisions. The hard 
parton level interactions are given by standard perturba- 
tive matrix elements and parton showers, which are not 
altered by the softer non-perturbative effects. The SCI 
model then applies an explicit mechanism where color- 
anticolor (corresponding to non-perturbative gluons) can 
be exchanged between the emerging partons and hadron 
remnants. The probability for such an exchange can- 
not be calculated and is therefore taken to be a constant 
given by a phenomenological parameter P. These color 
exchanges modify the color connections between the par- 
tons and thereby the color string-field topology, resulting 
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in different final states after the standard Lund model 
Jl2| has been applied for hadronization (Fig. p. 

The GAL model is similar in spirit, but is formu- 
lated in terms of interactions between the strings and 
not the partons. Soft color exchanges between strings 
also change the color topology resulting in another string 
configuration (Fig. Q). A generalization of the area law 
suppression e~ bA in the Lund model gives the probabil- 
ity for two strings to interact as P = Pq[1 — exp (—bAA)] 
depending on the resulting change AA of the areas swept 
out by the strings in momentum space. The exponential 
factor favors making "shorter" strings, e.g., events with 
gaps, whereas making 'longer' strings is suppressed. The 
fixed probability for soft color exchange in SCI is thus in 
GAL replaced by a dynamically varying one. 

The Monte Carlo implementations of SCI and GAL 
generate complete events with final state particles. This 
allows an experimental approach to classify events de- 
pending on the final state: e.g., gaps or no-gaps, lead- 
ing (anti)protons, charmonium etc. Thus, one obtains 
predictive models where a single parameter (P and Pq), 
regulating the amount of soft color exchanges, has a uni- 
versal value determined from HERA rapidity gap data. 

The SCI and GAL models give various diffractive hard 
scattering processes by simply choosing different hard 
scattering subprocesses in PYTHIA. Rapidity gap events 
containing a W, a dijet system or bottom quarks are 
found to be in agreement with Tevatron data The 
CDF dat a on dijets in DPE events are also repro- 
duced §, yfTboth in cross section and more exclusive 
quantities such as the dijet mass fraction. Thus, our 
models successfully pass these tests given by processes 
with similar dynamics as diffractive Higgs production. 

The properties of the Higgs boson in the Standard 
Model are fixed, except for its mass. The present lower 
limit is 114.1 GeV and \ 2 fits to electroweak data favors 
rn H < 212 GeV Q. The latest LEP data give an indica- 
tion (/~ 2.1 <j) of a Higgs with a mass of 115.6 GeV Jl6| . 
We therefore use run = 115 GeV as our main alternative, 
but also consider mn up to 200 GeV. 

Higgs production at the Tevatron and the LHC can 
proceed through many subprocesses, which are included 
in pythia version 6 111]]. The dominant one is gg — ► H, 
which accounts for 50% and 70% of the cross section (for 
115 < m H < 200 GeV) at the Tevatron and LHC, re- 
spectively. In this process, see Fig. |l|, the gluons couple 
to a quark loop with dominant contribution from top 
due to its large coupling to the Higgs. Other produc- 
tion channels are (focfe — > H, q^ — > Z H, qiqj — > W H, 
QiQj —* QkQlH and gg — > ququH- Their relative con- 
tributions depend both on the Higgs mass and the cms 
energy. The overall cross sections are obtained by fold- 
ing the subprocess cross sections with the parton density 
distributions (we use CTEQ5L ]Tij]). This basic factor- 
ization is proven for inclusive hard scattering processes. 
It is assumed to also hold in our model since the soft 
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FIG. 2: Predictions from the SCI model for Higgs produc- 
tion in single diffractive events denned by a leading proton 
or rapidity gap criterion at the Tevatron and LHC. Absolute 
cross sections and relative ratio (single diffractive to all Higgs) 
versus the Higgs mass. 



hadronization processes should not influence the cross 
section for the hard subprocess, but only affect the dis- 
tribution of hadrons in the final state. 

After the standard parton showers in PYTHIA, SCI or 
GAL is applied using the parameter values P=0.5 and 
Po=0.1, respectively. This gives a total sample of Higgs 
events, with varying hadronic final states. Single diffrac- 
tive (SD) Higgs events are selected using one of two cri- 
teria: (1) a leading (anti)proton with xp > 0.9 or (2) a 
rapidity gap in 2.4 < \rj\ < 5.9 as used by the CDF col- 
laboration. Applying the conditions in both hemispheres 
results in a sample of DPE Higgs events. The resulting 
cross sections and relative rates are shown in Fig. || as a 
function of the Higgs mass and for m# = 115 GeV in Ta- 
ble The results have an uncertainty of about a factor 
two related to details of the hadron remnant treatment 
and choice of parton density parameterization. 

The cross sections at the Tevatron are quite low in view 
of the luminosity to be achieved in Run II. Higgs in DPE 
events are far below an observable rate. For m# = 115 
GeV, only tens of single diffractive Higgs events are pre- 
dicted. Only the most abundant decay channel, H — ► bb, 
can then be of use and a very efficient 6-quark tagging 
and Higgs reconstruction is required. The conclusion for 
the Tevatron is that the advantage of a simplified recon- 
struction of the Higgs in the cleaner diffractive events is 
not really usable in practice due to a too small number 
of diffractive Higgs events being produced. 
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TABLE I: Cross sections at the Tevatron and LHC for Higgs 
in single diffractive (SD) and DPE events, using leading pro- 
ton or rapidity gap definitions, as well as relative rates (SD/all 
and DPE/SD) and number (#) of events, obtained from the 
soft color exchange models SCI and GAL. 
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In contrast, the high energy and luminosity available 
at the LHC facilitate a study of single diffractive Higgs 
production, where also the striking H — ► 77 decay should 
be observed. Also a few DPE Higgs events may be ob- 
served. The quality of a diffractive event changes, how- 
ever, at LHC energies. Besides the production of a hard 
subsystem and one or two leading protons, the energy is 
still enough for populating forward detector rapidity re- 
gions with particles. As seen in Fig. ||, the multiplicity of 
particles is considerably higher at the LHC, compared to 
the Tevatron. The requirement of a "clean" diffractive 
Higgs event with a large rapidity gap in an observable 
region cannot be achieved without paying the price of 
a lower cross section. Requiring gaps instead of leading 
protons gives a substantial reduction in the cross section, 
as seen in Table [| Note that the high luminosity mode of 
LHC cannot be used, since the resulting pile-up of events 
would destroy the rapidity gaps. 

The Monte Carlo model does not include any specific 
mechanism for the exclusive reaction pp — > ppH and our 
simulations did not produce any such events. 

For comparison we have also investigated single diffrac- 
tive Higgs production in the pomeron model. This is 
based on the Regge framework with the exchange of a 
Pomeron with vacuum quantum numbers p8[ , given by 
an effective pomeron flux [fij. In case of a hard scat- 
tering process, which resolves an underlying parton level 
process, a parton structure of the Pomeron may be con- 
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FIG. 3: Multiplicity (for LHC divided by 2.5) in the region 
2.4 < \r)\ < 5.9 in the hemisphere of a leading proton with the 
indicated minimum xf, for Higgs events from the SCI and 
the pomeron models. 



sidered |l[ and the data on diffractive deep inelastic scat- 
tering from HERA can be well described by fitting par- 
ton density functions in the Pomeron [|l9). Applying ex- 
actly the same model for pp gives, however, diffractive 
hard scattering cross sections that are up to two orders 
of magnitude larger than what is observed at the Teva- 
tron. Although this can be cured by appropriately mod- 
ified pomeron flux functions, it may indicate a deeper 
non-universality problem of the pomeron model || . 

To get numerical estimates we use the pomeron model 
implemented in the pompyt Monte Carlo ]20|] . The par- 
ton densities in the Pomeron are from a fit (parameteri- 
zation I in pl|) to the diffractive structure function mea- 
sured at HERA. The pomeron flux |22| has been renor- 



malized 23 so as to reproduce the observed relative rates 
of diffractive hard scattering processes both at HERA 
and the Tevatron. 

The pomeron model is constructed to give a leading 
proton with a spectrum essentially as 1/(1 — xp). It is 
developed for situations where xf — > 1 dominates and 
usually taken to be trustworthy only for xf > 0.9. As 
shown in Fig. |i|, however, this distribution is strongly 
distorted in this case due to the kinematical condition 
imposed by the Higgs mass. At the Tevatron energy, the 
cross section is dominated by smaller ip. This makes 
the results of the model sensitive to a phase space region 
where the pomeron model cannot be safely applied. In 
particular, the diffractive Higgs cross section will depend 
on whether the usual requirement xf > 0.9 is applied 
or not. The resulting cross section also depends strongly 
on what conditions for diffraction are applied. The re- 
quirement at the Tevatron experiments of no particles in 
the rapidity region 2.4 < |?7| < 5.9 imposes a very strong 
reduction. If the gap can be in a more forward rapid- 
ity region, based on extended detector coverage, a much 
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FIG. 4: Distribution in xf ~ p\\/pmax of leading protons 
in single diffractive Higgs production in the pomeron model 
(pompyt Monte Carlo) applied to Tevatron and LHC ener- 
gies. Curves are for all such events (dotted), events with no 
particles in 2.4 < < 5.9 (solid) and with a gap of size at 
least four units in rapidity (dashed). 



larger rate of diffractive Higgs is obtained as illustrated in 
Fig. [|. Similar, but not as strong effects are also present 
at LHC energies. 

In view of this, predictions for the diffractive Higgs 
cross section will be somewhat uncertain in the pomeron 
model. To give some numbers, nevertheless, we use cri- 
terion (1) with a leading proton with xp > 0.9, but no 
specific gap requirement. This gives a cross section of 
2.8 fb for single diffractive Higgs production at the Teva- 
tron and 410 fb at the LHC. This includes reduction fac- 
tors of 5.2 and 9.2, respectively, from the pomeron flux 
renormalization |2j| making HERA and Tevatron data 
compatible but leaving an extrapolation uncertainty for 
the LHC energy. 

In contrast to the pomeron model, the SCI and GAL 
models are constructed to describe different final states 
through a general mechanism for soft color exchanges 
giving a smooth transition between diffractive and non- 
diffractive events. This implies a better stability with 
respect to variations of the conditions used to define 
diffractive events. Moreover, the energy dependence of 
SCI and GAL has proven successful. Data on various 
diffractive hard scattering processes at HERA and Teva- 
tron are well reproduced. The soft color exchange models 
should, therefore, give more reliable predictions. 

In conclusion, we have investigated the prospects for 
discovering the Higgs boson in diffractive events having a 
lower hadronic background activity that should simplify 
the reconstruction of the Higgs from its decay products. 



We find that the rate of diffractive Higgs events at the 
Tevatron will be too low to be useful. Therefore, the 
Higgs must here be searched for in normal events with 
their larger hadronic activity. At LHC diffractive events 
are not as clean as expected, since the large available en- 
ergy produces an increased hadronic activity. Still, LHC 
should facilitate studies of Higgs in single diffraction and 
the observation of some DPE events with a Higgs boson. 



[1] G. Ingelman, RE. Schlein, Phys. Lett. B152, 256 (1985). 

[2] UA8 Collaboration, R. Bonino et al., Phys. Lett. B 211, 
239 (1988); A. Brandt et al, Phys. Lett. B 297, 417 
(1992); Phys. Lett. B 421, 395 (1998). 

[3] G. Ingelman, hep-ph/9912534. 

[4] M.G. Albrow, A. Rostovtsev, hep-ph/0009336. 

[5] A. Bialas, P.V. Landshoff, Phys. Lett. B 256, 540 
(1991); H.J. Lu and J. Milana, Phys. Rev. D 51, 6107 
(1995); D. Graudenz, G. Veneziano, Phys. Lett. B365, 
302 (1996); JR. Cudell and O.F. Hernandez, Nucl. 
Phys. B 471, 471 (1996); M. Heyssler, Z. Kunszt and 
W. J. Stirling, Phys. Lett. B 406, 95 (1997); V.A. Khoze, 
A.D. Martin and M.G. Ryskin, Eur. Phys. J. C 14, 525 

(2000) ; M. Boonekamp, R. Peschanski and C. Royon, 
Phys. Rev. Lett. 87, 251806 (2001); B. Cox, J. Forshaw 
and B. Heinemann, Phys. Lett. B 540, 263 (2002). 

[6] R. Enberg et al, Phys. Rev. D 64, 114015 (2001). 

[7] A. Edin et al, Phys. Lett. B 366, 371 (1996); Z. Phys. 

C75, 57 (1997). 
[8] J. Rathsman, Phys. Lett. B 452, 364 (1999). 
[9] A. Edin et al, Phys. Rev. D 56, 7317 (1997); C. Brenner 

Mariotto et al, Eur. Phys. J. C 23, 527 (2002) 
[10] G. Ingelman et al, Comput. Phys. Commun. 101, 108 

(1997). 

[11] T. Sjostrand et al, Comput. Phys. Commun. 135, 238 

(2001) . 

[12] B. Andersson et al, Phys. Rep. 97, 31 (1983). 

[13] CDF Collaboration, T. Affolder et al, Phys. Rev. Lett. 

85, 4215 (2000). 
[14] N. Timneanu et al, hep-ph/0206147, proc. DIS 2002. 
[15] D. Abbaneo et al, CERN-EP-2000-016, Jan 2000. 
[16] LEP Higgs Working Group, hep-ex/0107029. 
[17] H.L. Lai et al, Eur. Phys. J. C 12, 375 (2000). 
[18] K. Goulianos, Phys. Rept. 101, 169 (1983). 
[19] HI collaboration, T. Ahmed et al, Phys. Lett. B 348, 

681 (1995); C. Adloff et al, Z. Phys. C 76, 613 (1997). 

ZEUS collaboration, M. Derrick et al, Z. Phys. C 68, 

569 (1995). 

[20] P. Bruni et al, jittp : //www3 . tsl .mi. se/thep/pompyt/ 

[21] T. Gehrmann, W.J. Stirling, Z. Phys. C70, 89 (1996). 

[22] A. Donnachie, P.V. Landshoff, Phys. Lett. 191B, 309 
(1987); [Erratum-ibid. B 198, 590 (1987)]. 

[23] K. Goulianos, Phys. Lett. B 358, 379 (1995); K. Gou- 
lianos, J. Montanha, Phys. Rev. D 59, 114017 (1999). 



